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A mechanism using the position-dependent gauge coupling is proposed to localize non- 
Abelian gauge fields on domain walls in five-dimensional space-time. Low-energy effective 
theory posseses a massless vector field, and a mass gap. The four-dimensional gauge invari- 
ance is maintained intact. We obtain perturbatively the four-dimensional Coulomb law for 
static sources on the domain wall. BPS domain wall solutions with the localization mecha- 
5^ ■ nism are explicitly constructed in the (7(1) x 17(1) supersymmetric gauge theory coupling to 

CLi" the non-Abelian gauge fields only through the cubic prepotential, which is consistent with 

the general principle of supersymmetry in five-dimensional space-time. 
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§1. Introduction 



Besides supersymmetric theories,^ the most intriguing possibilities for unified 
theories beyond the standard model are models with extra dimensions, that is, the 
Q^' brane-world scenario.^ - ® In this scenario, our four-dimensional world is localized on 

defects (which are often called branes) such as domain walls in a higher-dimensional 
space-time, and the standard model particles are assumed to be localized on the 
defect. Domain walls are the simplest defect and have been most useful to construct 
realistic models. However, localization of gauge fields on domain walls has been 
notoriously difficult in field theories, although scalar and spinor fields have been 
successfully localized. It has been recognized early that the model with the warp 
factor does not help to localize the gauge field on the domain wallJ^® More recently 
it was proposed that bulk and boundary mass terms can be introduced into a warped 
O ■ model to help localize gauge fields!^ It was achieved at the cost of a subtle fine- 

tuning as well as certain boundary interactions to restore the gauge invariance. An 
explicit model of Abelian gauge field localized on a domain wall in five-dimensional 
space-time has been obtained using tensor multiplet.® However, it has been difficult 
to extend the idea of tensor multiplet to incorporate the non-Abelian local gauge 
symmetry. If we are content with toy models of lower-dimensional world-volume, 
such as domain walls with the three-dimensional world-volume, there have been 
a number of proposals, assuming nonperturbative effects,^' or using perturbative 
methodsP^'^ However, it has been difficult to obtain an explicit model of non- 
Abelian gauge fields localized on a domain wall in five-dimensional space-time. 

A basic problem has been pointed out to localize gauge fields on a domain 
wallP 1 ® We wish to obtain a (perturbatively) unbroken gauge symmetry on the 
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world volume of the domain wall. Since one wishes to prevent gauge fields to propa- 
gate freely in the bulk, one is tempted to consider that the bulk space-time outside 
of the domain wall to be in the Higgs phase where gauge symmetry is broken. Unfor- 
tunately, however, the flux coming out of the source on the domain wall is absorbed 
by the bulk in the Higgs phase and cannot reach beyond the width of the domain 
wall even in the direction along the world volume of the domain wall. Because of 
this screening effect, the vector field acquires a mass of the order of the inverse width 
of the wallP^'®®'^ On the contrary, if a vector field is confined in the bulk and 
deconfined on the domain wall, the flux coming out of a source should be expelled 
from the bulk, producing the four-dimensional Coulomb law on the world volume of 
the domain wall. Therefore we need to consider the confining phase for the bulk, 
whose explicit implementation is often difficult. This difficulty is particularly acute 
in our problem, since a realistic model requires the five-dimensional bulk, where the 
knowledge on nonperturbative effects is scarce. Nonperturbative effects in a five- 
dimensional gauge theory with a cut-off was proposed to obtain the layered phase 
that confines only along the extra dimension.^ 

Even if the detailed knowledge of nonperturbative effects is not available, the 
confining medium can be rephrased classically by introducing a dielectric permeabil- 
ity e for gauge fields EJiQSJ i n the classical electrodynamics of a dielectric medium, 
the electric flux density D is the sum of the electric field E and the polarization P 
induced by E 

D = e E + P = eE (14) 

where the permeability e(x) depends on the spatial position x, and the vacuum per- 
meability is denoted as eo- For ordinary dielectric media, the dielectric permeability 
e is greater than the vacuum permeability eo, since the polarization P is always 
induced in the same direction as the electric field E. On the other hand, it has 
been proposed that the confining vacuum can be represented by an unusual dielec- 
tric permeability e — > 0, namely by a perfect dia-electric medium. A relativistic 
version of the dielectric permeability can be expressed by a Lagrangian 14 ' with the 
field strength F^ u = d^W v — d u W^ of an Abelian vector field 

C = -\e{x)F, v F^. (1-2) 

The dielectric permeability in this form is nothing but the position-dependent gauge 
coupling g 2 {x) 

g 2 (x) 

and the region of the confining vacuum is represented by the strong coupling g 2 — > 
oo in this classical language. Therefore we can represent the confining vacuum in 
the bulk and deconfining vacuum on the domain wall classically by the position- 
dependent gauge coupling, by requiring strong coupling g 2 — > oo asymptotically in 
the bulk away from the domain wall, and weak coupling on the domain wall. 

The purpose of our paper is to propose the position-dependent gauge coupling 
as a model for a gauge field localization on domain walls, to examine its generic 
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properties, and to construct explicit examples of such domain walls as BPS states 
in supersymmetric gauge theories. The general features of the position-dependent 
gauge coupling are discussed in §2. We find that the four-dimensional gauge in- 
variance is intact, assuring the existence of massless gauge field in the low-energy 
effective theory. By analyzing modes of effective fields in four-dimensional world- 
volume, we find that our model indeed has a massless gauge field and a mass gap. 
It is worth stressing that this mechanism is applicable to localization of the non- 
Abelian gauge fields as well as the Abelian gauge field. We will discuss these points 
in §3. We also show in §4 that the gauge field localized on the domain wall exhibits 
perturbatively a four-dimensional Coulomb law for static sources localized on the 
world volume of the domain wall. As explicit examples, we construct BPS domain 
wall solutions in supersymmetric gauge theories with U(l) x U(l) gauge group in 
§5. The non- Abelian gauge fields couple to these U(l) x U(l) vector multiplets only 
through the cubic prepotential and are localized on the domain wall. It is remarkable 
that the cubic coupling among vector multiplets is just sufficient to give a nontriv- 
ial profile of position-dependent coupling function automatically, once the domain 
wall is formed as the background solution. This satisfies the stringent constraint of 
supersymmetric gauge theories in five-dimensional space-time allowing only up to 
cubic coupling in the prepotential of vector multiplets.^) We will discuss briefly 
mechanisms to introduce matter multiplets in the nontrivial representations of the 
localized non- Abelian gauge fields in §6. 



In this section, we will explore general features of localized gauge fields on do- 
main walls, without relying on a particular form of domain wall solutions. Concrete 
examples of domain wall solutions demonstrating the feasibility of our mechanism of 
localization will follow in later sections. We consider a domain wall in 4 + 1 dimen- 
sional space-time, whose coordinates are denoted by x M ,M = 0, 1, ••• ,4, whereas 
the world-volume coordinates and the codimension of the domain wall are denoted 
as x^, /i = 0, 1, 2, 3, and y, respectively. Namely the domain wall profile depends only 
on y. We assume that the five-dimensional gauge field acquires a position-dependent 
gauge coupling function e(y), on the background domain wall solution. We normalize 
the position-dependent gauge coupling function e(y) as 



and denote the four-dimensional gauge coupling as g^. Let us first consider gauge 
field Wt, with non- Abelian gauge group G, where o runs from 1 to dim G, in 4 + 
1-dimensional space-time. Denoting the source and the field strength as Jm and 
Fl IN = d M W^-d N W^-f abc W^W^ with M, N = 0, 1, ■ ■ ■ , 4, we take the following 
Lagrangian with the position-dependent gauge coupling function e(y) instead of the 
ordinary five-dimensional gauge coupling constant l/gf 



§2. A model of position-dependent gauge coupling 





C = --e(y)F MNa F\ 
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Our metric convention is T]mn = diag.(+, — , ■ • ■ , — ). 

We assume that the position-dependent gauge coupling function e(y) is real and 
nonnegative everywhere 

e(y) > 0. (2-3) 

and has a profile localized on a domain wall. Let us take the center of the domain 
wall to be at y = 0. We make a crucial assumption for e(y) to vanish at both infinity 

y — > ±oo 

e{y) -+ 0, y -> ±oo. (2-4) 

This asymptotic behavior represents the strong coupling in the bulk, namely the 
confining vacuum. It has been proposed that this type of "a perfect dia-electric 
medium" is the classical representation of confining vacuum El S3 
The Lagrangian ()2-2[) gives the field equation 

= e{y)(D N F NM ) a + (d»e(y))F$ M , (2-5) 

where (D N F NM ) a = d N F% M - f abc W Nb F^ M and we have used the fact that the 
position-dependent gauge coupling function depends only on y. By analyzing the 
field equation (|2-5|) . we will explore the low-energy effective field theory in four- 
dimensional space-time in the following sections. 

§3. Mode Analysis 

In this section, we would like to find massless and massive modes that appear 
in the low-energy effective field theory of the five-dimensional gauge theory with 
the position-dependent gauge coupling function e(y). We first need to obtain mode 
functions using the linearized field equation arising from quadratic terms of the 
Lagrangian (|2-2p . We can define a mode expansion of the Lagrangian assuming 
these mode functions are complete, provided they are normalizable. The linearized 
field equation is given by 

= e{y){d N d N W a M - d N d M Wfj) + (dye{y)){d y W a M - d M W%). (3-1) 

To find mode functions, we here choose an axial gauge 

W y a = 0. (3-2) 

Then the field equation (|3-1|) for M = \x component becomes 

= e(y) - d 2 y )W* - dpff'WS} - d y e(y)d y W«. (3-3) 

Operating by d^ 1 to the above field equation, we find 

= -d y {e(y)d y (d»Wp). 

Thus d^Wy separates variables x and y, and factorizes into 



QH W a = Y (y)F a (x), 



(3-4) 
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where we define 

Y(y)= f V dy'^-, (3-5) 



and F a {x) is a function of x only. Therefore we find no propagating modes in the 
longitudinal part d^W?. Even though the longitudinal part does not provide a 
propagating modes, it may contribute to a potential responding to a static source, 
as is usual in quantum electrodynamics. We will discuss the static potential due to 
a source in the next section. 

Now we decompose into transverse and longitudinal components, W^ T and 
W% . In particular, using (|3-4|) . the longitudinal component reduce to 

w aL _ J_q Qvy/a 



Y(y)±d,F a (x) (3-0) 



Plugging the decomposition = W" + IV" into (j3-3|) and using (|3-6p , we finally 
find the field equation for the transverse component W^ T 

= e{y)d v d v Wf - d y {e{y)d y Wf). (3-7) 

Let us now expand W^ T by a complete set of wave functions u n (y) in y as 

W^ T = ^2w a ^\x)u n (y), (3-8) 

n 

where w a ^ n \x) satisfies 

(d u d v + m 2 n )w< n \x) = Q, 

and m n represents mass of the mode. Plugging the mode expansion (j3-8j) into the 
field equation (|3-7j) . we find the equation for the mode function u n {y) 



1 d d 2 
e(y)—-m n 



e(y) dy dy 



My) = 0. (3-9) 



If we change the variable from y to Y defined in Eq. (|3-5p satisfying = e(y) J^, the 
mode equation f|3-9|) can be transformed into a bound state problem at the threshold 
(zero energy) for a potential U(Y) which is proportional to the mass squared of 
the mode 

H = -l-^ + U{Y) > U{Y) = - l -mle\y(Y)), (3-10) 
Hu n (Y) = 0. (3-11) 
Let us note that the normalization of the position-dependent coupling function e is 



fixed by the four-dimensional gauge coupling g^ in Eq. (|2-1|) 

i = J dye(y) = J dYe 2 (y(Y)), (3-12) 
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where the integration should be carried out over the entire region of y or Y. 

For a given shape e(y(Y)) of the potential, we can find all possible threshold 
bound states at various discrete depths of the potential U(Y) by adjusting rr? n . 
In this way, finding threshold (zero energy) bound state solutions gives a discrete 
spectrum of mass squared m 2 of modes. We wish to obtain the low-energy effective 
Lagrangian defined by an integration of the fundamental five-dimensional Lagrangian 
over y 

£ c ff = J dy£ 5 . (3-13) 

This effective Lagrangian dictates the normalization condition which naturally con- 
tains the following measure gl e(y) for the mode function 

dy [gl e(y)]u n {y)*ui{y) = 5 rd . (3-14) 

This measure is a distinctive feature of our threshold bound state mode functions. 
We see that constant mode is always a zero energy solution: 

uq = constant, Huq = 0. (3-15) 

The normalizability of this constant mode is equivalent to the condition of finiteness 
of the effective four-dimensional gauge coupling in Eq. (j3- 12|) . As will be illustrated 
by solvable examples in the following, there is a mass gap for these threshold bound 
states, and their wave functions are normalizable because of the nontrivial measure 
in Eq. (|3-14p . similarly to the constant mode, provided the finite four-dimensional 
gauge coupling can be defined by Eq. (|3-12p . This situation is quite different from 
the threshold bound states in the usual quantum mechanics problems in one spatial 
dimension. 

To illustrate our procedure of finding the spectrum of modes for the position- 
dependent gauge coupling, we will take two examples of solvable potentials that can 
also serve as approximations to our concrete examples of domain wall solutions in 
subsequent sections. 

Solvable Example 1 

We first consider the potential 

cosh 2 aY ' 

which is plotted in FigQJa). Because of the relation (|3T0p and the normalization 
condition (|3-12|) . the position-dependent coupling function e is fixed in this case as 



U(X) = ~ .._A (3-16) 



e(Y) = — ,p . (3-17) 

y J 54 V 2 coshaY v ; 

Then Eq. f)3- 5[) implies y = (g4\^2/a 3 ^ 2 ) sinh aY. Therefore, from Eqs. (|3-10p and 
(|3T5|) . we find 



£(y) 9 J 4coshaY(y) g A \j Ag\ + 2g A a^ ' (348) 
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which is shown in FigQTb). 

If we consider the eigenvalue problem Hvk = E^Vk, we immediately fincP^ 
that the finite number of bound states exist with a discrete energy spectrum E)~ = 
-{a 2 /2)(s - k) 2 , k = 0, 1, . . . < s with 




s = - -1 + A /! + -£ (3' 19 ) 



The threshold bound state E^ = occurs if and only if k = s with s being a 

nonnegative integer s = n, n = 0, 1, 2, This gives the n-th threshold bound state. 

In that case, the potential depth Uq satisfies 

a 2 

U = —n{n + l). (3-20) 

Eq. (|3T0p implies the threshold bound state spectrurrO 

ml = 2gjan(n + 1), n = 0,1,2,.... (3-21) 

First we consider the energy level at n = 0. This gives the massless mode m\ = 0, 
whose wave function turns out to be the constant uq(u) = 1, as we have seen before. 
Secondly, the first excited mode gives m 2 = \g\a. Similarly higher excited modes 
gives larger discrete values of m\. The normalizability (|3T4p of the n-th mode in 
this example is given by the finiteness of 

dye{y)\u n (y)\ 2 = / dY — -5— K(Y)| 2 . (3-22) 

) J -00 2g| cosh aY 



*' The continuum spectra with positive energy Ek do not contribute, except the limiting case of 
zero energy Ek — > 0. By regularizing in a finite interval in Y", we find that the zero energy solution 
reduces to our n — solution in the limit of infinite interval. 
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We find that wave functions of all the excited modes are normalizable, since they are 
just polynomials^ in tanh aY. Thus we conclude that there is a massless mode and 
a mass gap for the first excited mode, both of which are normalizable. We can safely 
use the effective field theory of the massless gauge fields below the energy scale mi, 
ignoring the massive modes. The value of the mass gap is proportional to the inverse 
of the width ^ of the position-dependent coupling function e and the square of the 
gauge coupling g\ . The massless zero mode wave function including the square root 
of the measure y e(y)uo(y) is plotted in FigQJb) by a dashed line. 

Solvable Example 2 

Next solvable example is a square well potential 

U ( Y 1 = { o U °' |y| > 1 ■ < 323 » 

(See Figj2j) Then, the coupling function is given by 



/2U0 



111 < a 

ilhfl ■ I I 

0, \Y\ > a 



Using the normalization of the profile e(y), we determine 

Energy levels and wave functions for this square well potential can be exactly solved, 
and the threshold bound states occur when we choose 



^o = ^ 2 , n = 0,1,2,.... (3-24) 
8a z 



Therefore the mass spectra of threshold bound states are given by 



TT 



2„2 



m n = ^94,aU = —g 4 n 
la 

Again we find there is a massless zero mode at the level n = and finite mass gaps 
for higher excited states. 

§4. Four-dimensional Coulomb law 

In this section, we wish to demonstrate that the position-dependent gauge cou- 
pling in Eq. (j2-2p exhibits the four-dimensional Coulomb law for the static source on 
the world volume of the domain wall, following a treatment in Ref JIHp . We introduce 
a source term to Eq. t|2-2|) to examine the response 

£ = -\e(y)F a MN F aMN + J a M W aM . (4-1) 
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Fig. 2. The square well potential with a depth Uo- 
choose Uo suitably. 



There exists a threshold bound state if we 



We are interested in an external source Jm(x, y) that is localized on the world volume 
of the domain wall, and has components only in four-dimensional world-volume: 

JSt(x,y) = S(y)S M ^(x). (4-2) 

The free field equation for the gauge field W^, namely the equation of motion 
ignoring the nonlinear interaction terms reads 

= e(y)(d M d M W% - d M d N Wfo) + dye{y){d y W a N - d N W%) - J a N . (4-3) 

In this section we choose the Lorentz gauge (in five dimensions) 

8 M W^ = 0. (4-4) 

Then the field equation becomes 

= e{y)d M d M W% + d«e(y)(ayW% - d N W£) - J%. (4-5) 

Since the source has no y-component, the field equation for the extra-dimensional 
component W y becomes 

= e(y)d M d M W£, (4-6) 

without any source term at y = 0. Because of nonnegative definiteness of e(y), we find 
that Wy obeys a free field equation in five dimensions without source. Assuming that 
there is no external source at infinity, we obtain that there is no nontrivial solution. 
To demonstrate it, let us solve the free field equation by going to the Euclidean space. 
Denoting the mixed representation of the momentum space only in four dimensions 
as Wy(p,y), we obtain 

= (p 2 -d 2 y)W y i (p,y). (4-7) 
There are two independent solutions with C+ (p) , C_ (p) as arbitrary functions of p 

W$(p, y) = C + (p)ePy + C_(p)e-™, (4-8) 
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which is valid in the entire region — oo < y < oo. Since no external source at infinity 
requires Wy(p,y) — > at both infinities y — > ±00, we obtain C + (p) = C_(p) = 0, 
which implies no nontrivial solution: Wy = 0. 

Taking account of Wy = and Eq. fj4-2|) , we obtain the field equation for W" as 

= e{y)d v d v W« + du {e(y)d y WZ) - 6(y)J^x). (4-9) 

Now we give a simple example to see the Coulomb law on the domain wall. We 
assume the weak coupling region (domain wall) is sufficiently thin and the coupling 



_ %! 

PC 

totic bulk gauge coupling 



profile behaves as e t hin(y) 

94 

For regularization purposes, we will add a large, but finite values for the asymp- 



ewn(y) = 6 -^ + \. (4-10) 

at at 

We consider this simplified situation of the zero-width limit of the domain wall 
in Eq. (|4-lU|) . and examine the case of the finite width of the domain wall later to 
confirm that our result is unchanged. Going again to the Euclidean space, and using 
the mixed representation W^(p,y) of the momentum space only in four dimensions, 
we obtain the field equation as 



( ^2_ d 2 ) 5(y)p 2 W-(p,y) dy(5{y)d y W*(p,y) 

= T J ^W^p,y) + f * 2 ^-%)J^(p). (4-11) 

#5 34 3i 

Since the field equation (|4-lip for y 7^ is identical to the Wy component, we obtain 

K(p,y) = C« + (p)e-Py0(y) + Of {p)^9(-y). (4-12) 

Integrating the field equation (|4-lip in the infinitesimal interval between y = —e and 
y = e, we obtain 

f -e d 2 W a (v v) n 2 ~ 1 f £ / ~ \ 

= _ / dy y f^lMl + *W2(p,v) - - / dyd y (5(y)d y W»(p,y)) - J%{p). 

(4-13) 

By inserting the solution (|4-12p for y 7^ 0, we find that the absence of divergent 
terms such as (S(y)) 2 in Eq. (|4-13p from the third term S(y)d y W®(p,y) requires 

C a +{p) = C a -(p) = C*(p). (4-14) 

With this condition (|4-14p . the third term vanishes^ : 

- 1 f dydy (5{y)dyW«{ P ,y)) = f dyd y (s(y)d y (sgn(y) e -^)) = 0. 

Q A_ 1/ — £ 9 A^ J — S 

(4-15) 

*' For an arbitrary smooth function f(y), we find dy f (8(y)sgn(y)) = 

- d y^ 5 (y) s S n (y) = °» where sgn(y) = y/Jy\. One can demonstrate it by using 5(y) = ^sgn(y) 
and a regularization such as the e — > limit of sgn(y) = tanh(y/e), or sgn(y) = y/(2e) for — e < y < e 
and sgn(i/) = sgn(j/) otherwise. 



Non-Abelian Gauge Field Localized on Walls 



11 



Now the field equation determines the amount of discontinuity 2pC^ (p) of the deriva- 
tive of W® at y = in terms of the source J„(p) leading to 

K(p,y = 0) = 9 \ J». (4-16) 
V 2 + -tP 

In the strong coupling limit g 2 — > oo of the bulk asymptotic coupling, we finally 
obtain 

km W^(p,y = 0) = §J^p). (4-17) 



If we put a static charge as the source J^{p) = Q a 5^,o, we obtain the potential in the 
coordinate space by a Fourier transformation of (|4T7p 



W«(x,y = 0) = ~g 2 Q a 5^ (4-18) 

with the three-dimensional spatial distance r. Thus we obtain the Coulomb law in 
four-dimensional world volume as we anticipated. 

The intermediate form of our potential in Eq. (|4-16p turns out to be identical to 
the result in Ref ll8p . However, let us note two important difference of our analysis 
from that in Ref.18). Firstly, we have introduced the bulk asymptotic coupling 
g 2 merely as a regularization parameter, and the agreement of the potential at an 
intermediate step is somewhat technical. Moreover our starting Lagrangian possesses 
the five-dimensional field strengths F^ N localized at the wall, in contrast to their 
Lagrangian in Ref, 18) where only four-dimensional field strengths F® are localized. 
This difference results in the presence of the third term in Eq. (|4-lip . As we noted, 
this term gives us a consistency condition (|4-14p . which does not follow from the 
field equation in Ref ll8p . Instead they seem to have assumed (quite naturally) a 
symmetry of W%(p, y) under y — > —y. Thanks to this consistency condition, the 
third term in Eq. (|4-13|) vanishes and the resultant potential has become identical to 
that in RefHSJ. 

So far we have been studying the zero-width approximation for the domain wall. 
Let us now examine the case of the finite width of the domain wall by regularizing 
the delta function profile. As a simplest regularization, we take the following step- 
function ansatz^ for the domain wall profile function e(y) (see Figj3]): 

Since the solution for Wy in Eq. (|4-8p uses only the positive definiteness of the cou- 
pling without referring to the y-profile of the coupling in Eq. (|4-6p . we find that 
Wy = even for the finite width case. Moreover, the position-dependent gauge cou- 
pling factors out for y ^ 0, and the source term exists only at y = 0. Therefore the 
solution for y ^ as well as the discontinuity at y = are unchanged from the zero- 
width case. We thus find that exactly the same solution (|4-17|) is valid in this finite 
width case. Generally we should obtain somewhat different y-profile of the solution 
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Fig. 3. The position-dependent coupling is very thin step function. In the limit of e — > 0, the 
position-dependent coupling becomes a delta function. An asymptotic height of l/gi is intro- 
duced for a regularization. 



for other finite width regularizations. However, our example shows that the quali- 
tative behavior, the four-dimensional Coulomb law particularly, should be the same 
as the zero- width case. In the case of Ref J8j) . the finite width regularization by the 
step function gives a solution for W% different from the zero- width limit, contrary 
to our Lagrangian. This result arises from the fact that our Lagrangian contains 
the extra-dimensional component of the field strength with the position-dependent 
coupling e{y). 

§5. Supersymmetric models for BPS walls 

5.1. General set up 

In order to have a realistic brane world with four-dimensional world-volume, we 
construct a domain wall using five-dimensional supersymmetric gauge theories, which 
have eight supercharges and consist of vector multiplets and hypermultiplets. We will 
consider at least two U(l) vector multiplets labeled by I = 1, 2, • • • , which contain 
gauge fields W^, M = 0, 1, • • • ,4 and neutral scalar fields U 1 , besides fermions. The 
vector multiplets for the non-Abelian group G with the dimension dim G also contain 
gauge fields Wjfa and scalar fields E a in the adjoint representation a = 1, • • ■ , dim G. 
Hypermultiplets as matter fields contain scalar fields Ha (besides fermions) with A 
labeling different flavors of hypermultiplets. Non-vanishing values of these Ha will 
break the U (1) gauge symmetries and give domain wall solutions. Since we do not 
wish for non-Abelian gauge symmetry to be broken, matter scalars H a are assumed 
to be singlets of the non-Abelian gauge group. 

It is easy to construct domain walls with a number of hypermultiplets inter- 
acting with vector multiplets, provided the gauge group involves one or more U(l) 
factors allowing the Fayet-Iliopoulos (FI) termP^'^^^ For simplicity, we use the 
U(l) gauge groups to build domain wall solutions, and the minimal kinetic term for 
U(l) vector multiplets. We choose all the FI terms along the common direction in 
the SU(2)r space. We will also use the strong coupling limit of these U(l) gauge 
couplings, whenever we wish to give an explicit exact solution of domain walls. 



Non-Abelian Gauge Field Localized on Walls 



13 



5.2. Domain wall sector 

We consider Np charged hypermultiplet^)] Ha, A = 1, • • • , Np with the charge 
qf,I = 1,2 for the U(l)\ x £7(1)2 gauge group. Neglecting the non-Abelian gauge 
group in this subsection, we obtain the bosonic part of the Lagrangian for the domain 
wall sector 

£wall = ~ {Fl IN f + ^(OmU 1 ) 2 + \V M H A \ 2 - V, (5-1) 

where T>mHa = (Om + iWf^qf)H A is the covariant derivative. The potential V is 
given by 

V = KqfZ 1 - m A )H A \ 2 + ^{Y 1 ?, (5-2) 

where rriA £ R is a real mass of the A-th hypermultiplet. The auxiliary fields Y 1 are 
given by their equations of motion in this model as 

Y 1 = ej (cj - qf\H A \ 2 ) , (5-3) 



where the FI parameteJ**^ for J-th U(l) factor group is denoted as cj. 

Taking a Bogomolnyi completion, we obtain the Bogomolnyi bound which is sat- 
urated by the Bogomolnyi-Prasad-Sommerfield (BPS) equation. The energy density 
is given by 



£=-L(d y E I -ej( Cl -qf\H A \ 2 )) 2 



2ej 



+ \V y H A + (qfZ I -m A )H A \ 2 

+8 y (aZ 1 - qfZ^HAl 2 + m A \H A \ 2 ) , (5-4) 

where we have chosen the gauge Wy = 0. The topological charge for the domain 
wall connecting the vacuum /3 to a can be read from the last term to give the tension 

/+0O 
dydyf, with / = aU 1 - qfz'lHAf + mAlHA] 2 . (5-5) 
-oo 



The BPS equations are obtained from Eq. fj5-4[) ai 



(dy + qfiS 1 + iWy r ))HA = H A m A , (5-6) 
dyS 1 = ej( Cl -qf\H A \ 2 ). (5-7) 



*' Although a hypermultiplet contains two complex scalars H\,i = 1,2 for each flavor A, we 
denote only one of them H\ as Ha, since the other one H A does not participate in our BPS domain 
wall solutions (H A = 0). 

**' Both the FI parameter c aI and the auxiliary fields Y aI , a = 1,2,3 are triplets of SU(2)r, 
as in Eq. (|5-2I|) . We have chosen the direction of FI parameters for the both U(l)i, I = 1, 2 to be 
parallel along the third direction c aI = (0, 0, ci) and suppress to write components of auxiliary fields 
except the third component which we denote as Y aI — (0, 0, Y 1 ). 
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The gauge group indices / are summed for a fixed flavor index A in the first BPS 
equation (|5-6p . and vice versa in the second BPS equation (|5-7p . When these equa- 
tions are satisfied, the energy density becomes equal to [f]y=J^- 

The first BPS equation (|5-6p can be solved in terms of a constant matrix called 
moduli matrix Hqa (in our case of U (1) gauge theory, Hqa is actually a vector) EHUD 

H A = H QA Q^ l2 Q~ q ^ l2 e mAV , (5-8) 

where Qj,I = 1,2 are given by the solution to the master equation BZM22J which gives 
the solution to other BPS equation (|5-7p . In the strong coupling limit e\ — > oo, we 
can find exact and explicit solution fij from the following algebraic equations^ 

CJ = \H 0A \ 2 Q~ q * Q- q *e 2mAV . (5-9) 

In terms of these f2j, the hypermultiplet scalars are obtained by Eq, (|5-8p . whereas 
the vector multiplet scalars S 1 are given by 

U 1 = ^OylOgQj. (5-10) 

We first consider the case of four hypermultiplets which allows an easy construc- 
tion of appropriate domain walls, and then the case of three hypermultiplets as a 
model with the minimal number of matter fields. 

5.2.1. Two copies of two charged matter fields (four flavor model) 

The simplest model with a domain wall is the U(l) gauge theory containing two 
charged hypermutiplets with different masses We will just take two copies of such 
models. The charges qf, 7 = 1,2 for gauge group U(l)i and masses tua of the ^4-th 
hypermultiplets are given in TablelH We also choose FI parameters to be positive 





H A =i 


H A=2 


H A =3 






for C/(l)i 


1 


1 










for U{1) 2 








1 


1 




m A 


m 
2 


m 
— 2 


m 
2 


_ 2 



Table I. Charge and mass of matter fields (hypermultiplets) of the four-flavor model 



Ci = C2 = C > 0. 

The BPS domain wall solution is well-known for the two flavor model. We take 
the strong coupling limit ef — > oo where the model reduces to a nonlinear sigma 
model with the T*CP l target space, allowing an explicit exact solution for a domain 
wall from Eq. (j5-9p . By choosing the boundary condition for U —> =pm/2 at y — > =Foo, 
we find 

e f(y-yi) e =^{y-yi) 
Hi = a/c ==, H 2 = yfc- 



A/2coshm(y — y\) ^2 cosh m(y — y\) 

£ = — tanhm(y - yi), (5-11) 
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where the physical meaning of the moduli parameter y\ is the domain wall position. 
Precisely the same form of solution is valid for the second copy of the £7(1) model, 
with a moduli y 2 for the position of another wall 



U 2 



e 2 



(y-V2) 



m 



^2 cosh m(y - y 2 ) 
tanhm(y — y 2 )- 



H 4 = yfc 



e 2 



K3/-3/2) 



^2 cosh m(y - y 2 ) 



(5-12) 



If we take a difference of these two neutral scalars E 1 , we find a profile suitable 
for the position-dependent coupling function e(y), provided y\ < y 2 



U 1 -S z 



jrt 

— (tanhm(y 



y x ) - tanhm(y - y 2 )) , 



(5-13) 



which is positive definite, and falls off exponentially fast towards both infinities 
similar to the profile of the position-dependent coupling discussed in §3. 

The biggest advantage of this model is its simplicity. The profile of U 1 — U 2 is 
positive definite for y 2 > y\ and has a three layer structure with two outer skin with 
the width 1 /m and inner wall with the width y 2 — yi, as shown in FigUJ Therefore we 
can choose arbitrary wall width by adjusting the moduli y 2 ,yi whereas the domain 
wall skin width is fixed by the mass parameter of the model m. The wall profile can 
be as close as the step function, by choosing the mass parameter large m — >• 00, with 
a fixed y 2 — y\. On the other hand, the model becomes unstable for y 2 < y\ because 
of negative kinetic term for gauge fields. For that reason, one may be tempted to 
consider an another model with no moduli for the adjustable wall width, to which 
we turn next. 

5.2.2. U(l) x U(l) model with three Higgs flavors 

As is most easily seen by taking the strong coupling limit (|5-9|) . each gauge 
group acts as a constraint on hypermultiplet scalars to form domain wall solutions. 
Therefore the minimal number of flavors to have domain wall solution is the case of 
3 charged hypermultiplets Ha, A = 1,2,3. 

Let us take a model containing hypermultiplets with the charges qf, I = 1.2 and 
masses itla as given in Table ITU We easily find two supersymmetric vacua* a The 





H A =i 


H A =2 


H A =3 


U(l)i 


1 


1 





[7(1)2 





-1 


1 


VtlA 


m 









Table II. Charges and masses of the three matter fields 

first vacuum is given by 

#1 = 0, H 2 = ^c\-, H 3 = ^cTT^, Z 1 



0. 



(5-14) 



This model has another supersymmetric vacuum, which will not be used in this paper: 
the vacuum with Hi = y/ci, H3 = 0, S 1 — E 2 = m. For the second matter field, the second 
complex scalar of the hypermultiplet H^ 2 has a vacuum value H^ 2 = ^fSi instead of the first one 
H2 = HiT 1 = 0, similarly to the model considered in Refl27p. 
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E 1 , E2 


'E 1 - 


E 2 
































y 


y 









(c) (d) 

Fig. 4. A plot of the profile of the domain wall with two Higgs flavors. The shape of the difference 
E 1 — S 2 is similar to the position-dependent coupling discussed in §3 if \yi — J/2I ~ 1/m. ((a) 
and (b)) The shape also tends to the step-function if \yi — 2/2 1 2> 1/m. ((c) and (d)) 



The second vacuum is given by 

H x = ^cT, H 2 = 0, H 3 = V^, Z 1 = m, U 2 = 0. (5-15) 

Without loss of generality, we can choose the moduli matrix to be Ha = 
(y/ci, y/c±e myo , y/c%). The moduli parameter yo is taken to be real Ell- The BPS 
domain wall solution connecting these two vacua is given in terms of the solution of 
the master equations Q\ , Q 2 as 

Hi = yfce my f2- 1/2 , H 2 = ^r l e m y a n- l/2 n\ 12 , h 3 = ^q~ 1/2 , (5-ie) 

where Q\ , Q 2 are given in the strong coupling limit as 

Q l = e 2m y + e 2myo Q 2 , (5-17) 

1 _ e 2m(y-y ) + I _ e 2m(y-y ))2 _|_ 4(1 + £l) e 2m(y-y ) 

th = V 2(1 + a) ~ ■ (5-18) 

V C2 ' 

The vector multiplet scalars 17 1 , E 2 are given by (|5- 10)) . Since the domain wall 
solution connects two vacua in (|5-14p and (|5-15p . U 2 is appropriate to give the 



*' The generic moduli yo is complex, whose real and imaginary parts correspond to the domain 
wall position and the relative phase of two vacua. Since we are not interested in the phase moduli, 
we take j/o to be real here. 
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> 


's 2 


E 2 




y 


y 



(a) (b) 



Fig. 5. A plot of the profile S 2 of the domain wall with three Higgs flavors. The height of the 
function is suitably normalized to compare the profiles of different FI parameters with each 
other, (a) The shape is similar to the position-dependent coupling in §3 when C1/C2 ~ 1. (b) 
We can also obtain the step- function like profile when C1/C2 3> 1. 



position-dependent coupling function in this model. We easily find that the profile 
of E 2 is reflection symmetric with respect to the domain wall position y = y$. Let 
us take the domain wall position at the origin yo = 0. We obtain 

fh(-y) = (^(y)) -1 , ^ 2 (-y) = z\v)- (5-19) 

The asymptotic behavior at large values of y is given by 

U 2 (y) ^me- m{y ~^\ Ay = l ° g ^\ ( 5 . 2 0) 

m 

We see that the domain wall has a three layer structure: the overall wall width Ay is 
proportional to log(ci /C2), whereas the outer skin has the width 1/m. In this model, 
both the outer wall (skin) width and the inner wall width are fixed by parameters 
of the model, and are not the moduli of the domain wall. The exact profile given in 
Eqs. (|5-18p and (|5T0p is illustrated for different values of c\/c2 in Figj5j 

5.3. Position- dependent coupling function from the cubic prepotential 

Interactions between hypermultiplets and vector multiplets are specified by charge 
assignments of the hypermultiplets, whereas interactions among vector multiplets are 
specified by the so-called prepotential a(U). It has been found from general princi- 
ples® that the prepotential, which gives the Chern-Simons coupling for gauge fields 
together with other terms, in our five-dimensional theory should be at most cubic in 
vector multiplets. Let us write the bosonic part of a Lagrangian with the prepoten- 
tial a(U), by denoting the group label and gauge generators collectively as /. Each 
U(l) factor group can have a triplet of the FI term with the parameters c aI and the 
auxiliary fields Y aI , a = 1,2,3. Restoring two complex scalar H irA ,i = 1,2 for the 
hypermultiplets with the color (G) indices r and the flavor indices A, we obtain 

C = aiJ (-lF[f N F JMN + ^D M E I D M E J + l -Y aI Y aJ ^ - c aI Y aI 

+ aiJK [~y/ MNPQw I (fmnFpq + \[W M , W N ] J F* Q + ±[W M , W N ] J [W P , W Q ] K ^j } 
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+ (V M H irA )*V M H irA - (H^f^E 1 - m A f] r s H isA 

+ (H irA )*(o a ) i j q I (Y aI ) r s Hi aA , (5-21) 

where tua is the mass of the A-th hypermultiplet, and the derivative of prepotential 
a(£) is denoted by subscripts like 

_da(E) _ d 2 a(E) _ d 3 a(Z) 99 , 

Covariant derivatives T>m and q 1 U 1 — itia are understood to contain both U(l) 
and non-Abelian components with appropriate charges or representation matrices. 
We note that the gauge field kinetic term multiplied by the scalar E 1 arises as a 
supersymmetric completion of the Chern-Simons term, both of which follow from 
the cubic prepotential. 

The minimal kinetic term for U(l) vector multiplets is represented by a term of 
the form U 1 U 1 in the prepotential. As given in the previous subsections, the domain 
wall solution leads to a nontrivial kink profile for hypermultiplet scalars H A , and the 
vector multiplet scalars E 1 . Let us call those hypermultiplets and vector multiplets 
participating to form the domain wall as the domain wall sector. As described in 
§1, we wish to avoid the bulk in the Higgs phase in order to obtain localized gauge 
fields. This requirement is achieved by demanding the domain wall hypermultiplets 
to be neutral under the non-Abelian gauge fields which we wish to localize. Then 
the non-Abelian gauge fields can couple to the domain wall sector only through the 
prepotential among vector multiplets. Since the non-Abelian vector multiplets can- 
not appear linearly, the coupling between non-Abelian vector multiplets and domain 
wall sector should be linear in (a linear combination of) U(l) vector multiplets in 
the domain wall sector. By choosing a linear combination a\E l + a^E 2 of two ?7(1) 
vector multiplets, we can obtain a desired profile of the position-dependent cou- 
pling function e{y) at both infinities y — > ±oo, namely the asymptotically vanishing 
profile with a peak in the middle. Therefore we shall consider gauge group to be 
U(l) x U(l) x G with G as the non-Abelian gauge group which we wish to localize 
on the domain wall, and assume the following prepotential 



a (S) = ^(r 1 ) 2 + ^(i7 2 ) 2 + ^(mS 1 + a 2 U 2 )U a U a (5-23) 



whose constant coefficients a\ , a 2 are chosen appropriately for each model of the 
domain wall sector. The first and second terms reproduce the minimal kinetic terms 
for the U{1) vector multiplets, as we assumed in previous subsections. 

Since the non-Abelian vector multiplet U a appear only quadratically in the 
prepotential because of gauge invariance, it is easy to see that the above prepotential 
allows the BPS domain wall solution in the previous sections to remain a solution to 
the entire system of field equations. Therefore we can safely choose the BPS domain 
wall solution as the background solution and consider the effective Lagrangian on 
the domain wall. If we choose the following coefficients of the prepotential, we obtain 
the position-dependent gauge coupling function e{y) for the non-Abelian gauge fields. 
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For the model with four matter fields in §4.2.1, we choose 



di = — Q2 = a > 0. 



(5-24) 



For the model with three matter fields in §4.2.2, we choose 



oi = 



a 2 > 0. 



(5-25) 



The value of the effective gauge coupling in four-dimensional world- volume in Eq. (|2-ip 
can be adjusted by choosing the value of these coefficients a or 02. 



In this paper, we have discussed the localization of the massless vector fields by 
means of the position-dependent gauge coupling. We gave a few concrete examples 
where there exist the massless vector field and mass gap, and the Coulomb law 
emerges inside the domain wall. We expect that these localization properties do not 
depend on the details of the concrete coupling functions. If the position-dependent 
coupling function is everywhere nonnegative and vanishes at both infinities (weak 
coupling only at the center of the domain wall), the localization of the gauge field 
should occur in a similar way. 

The position-dependent coupling function desired for the localization can be 
realized by the cubic prepotential of the five dimensional supersymmetric gauge the- 
ory. The coupling function comes about thanks to the profile of the domain walls 
of the Abelian subsectors. We expect that these coupling functions offer appropri- 
ate examples for the localization of non- Abelian gauge fields, although the explicit 
functional forms of the coupling functions in our concrete examples of domain wall 
are somewhat more involved than our solvable examples in §3. 

In this paper, we considered the tree-level prepotential, in order to obtain the 
position-dependent coupling. The prepotential of the supersymmetric gauge theory 
with eight supercharges receives the nonperturbative quantum corrections generally. 
It is an interesting future problem to explore if our mechanism of localization of 
gauge fields due to the position-dependent coupling may be realized as a result of 
the nonperturbative quantum effects. 

We have succeeded to localize non- Abelian gauge fields. However, we still need to 
introduce matter fields in nontrivial representations of the non- Abelian gauge group, 
in order to build the standard model localized on the domain wall. It has been found 
that a non-Abelian flavor symmetry for degenerate hypermultiplets provides non- 
Abelian orientational moduli for domain walls. 30 ' These orientational moduli arises 
in nontrivial representations of the non-Abelian flavor group. If we promote (a part 
of) the flavor symmetry to a local gauge symmetry, these non-Abelian orientational 
moduli fields become matter fields interacting nontrivially with the non-Abelian 
gauge fields. We can introduce matter fields coupled to the non-Abelian gauge fields 
of our model in this way. In order to construct the standard model localized on the 
domain wall, it remains to see if matter fields in appropriate (chiral) representations 
can be introduced into our framework. The interaction of localized matter fields 
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and the localized gauge fields offers an intriguing question of charge universality.^ 1 
Our model should provide a concrete example of localized matter fields assuring the 
charge universality. 

We have succeeded to embed the BPS domain wall solutions in flat space into 
the five-dimensional supergravity theory, and found a model with the warped extra 
dimension j s a \ so an interesting future problem to embed our mechanism of 
localized gauge fields into supergravity. In this context, it is worthwhile to examine 
a recent objection against the Fayet-Iliopoulos parameter in super gravity.^ 
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Appendix A 

Threshold bound states and wave functions 



In this Appendix, we consider an eigenvalue problem of the following static 
Schrodinger equation discussed in §3 



1 d 2 



2dY 2 



+ U(Y) 



i/)(Y) = Eip(Y) for 



oo < Y < oo. 



(Al) 



We will find wave functions of the threshold bound states satisfying E = by 
choosing a suitable height of the potential U(Y). 

A.l. cosh potential 

Let us first take the potential U(Y) = —Uq/ cosh 2 aY in Eq. (|3T6|) (see FigJT]). 
Changing the variable by £ = tanhal" (—1 < £ < 1) and defining 2Uo/a 2 = s(s + l), 
and £ = \/—2E/a, (|A-1[) becomes 



d£ 



■i-e 



+ s(s + l) 



m = o. 



(A-2) 



In order to find a series solution around £ = 1, we define u = (1 — £)/2 and assume 
if; = u 1 (f>(u) = YlbLo a l u ■ By requiring a regular solution at £ = 1, we find 7 = £/2. 
Using <f>(u) instead of tft, Eq. (|A-2| ) reduces to 



d 2 d 
u(l - u)^2 + 1 + 5(1 - 2u)— - (S - s)(£ + s + 1) 



<t>(u) = 0, 



(A3) 



where <j){u) = YlhLo a i ul - ^ n order to obtain solutions regular at u = 1 (£ = —1), the 
series in 4>(u) = ^ aiu 1 has to terminate and £ = s — n. Thus the n-th eigenfunction 
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is given by the following n-th order polynomial 4> n (u) with < n < s for a given s 

, . . ^ 1 (l + 2s-n)\ n\ (s-n)\ , 
0n(n) = a ^ - {2s _ n)[ {n _ l)l{l + s _ n)l (- u ^ ■ 



1=0 

Threshold bound state is given by £ = 0, that is, s = n, which is achieved by 
choosing the potential height Uq = a 2 n(n + l)/2 as in Eq. (|3-20p . The wave function 
of the n-th threshold bound state is given by 

Mu) = aojr-^^±^(- U y. (A4) 
^ {i\y (n-l)l 

A. 2. square well potential 

Here we consider the square well potential with a finite potential step in Eq. (|3-23p 
(see FigJ2J. The reflection symmetry Y —> —Y dictates that eigenfunctions must be 
either an even or odd function. For \Y\ < a, we obtain the eigenfunctions ip(Y) = 
A cos kY for even functions, and B sin kY for odd functions, with k = \J2(E + Uq). 
For Y > a, the boundary condition at infinity requires that tp(Y) = Ce-^ Y \ with 
(3 = \/—2E. The bound state exists if and only if E < 0. The connection condition 
at Y = a gives k tan ka = f3 for even functions and k cot ka = —f3 for odd functions. 

The threshold bound state with E = means that j3 = 0, that is, the wave 
function becomes a constant outside of the well. The connection condition for /3 = 
is solved by discrete values of A:: k n = §^ with nonnegative integer n, where n is 
even for the even wave function and odd for the odd wave function, respectively. We 
obtain the n-th threshold bound state by choosing the potential depth Uq = k 2 /2 = 
n 2 7r 2 /(8a 2 ) as in Eq. (|3-24p , The wave function of the n-th threshold bound state is 
given for even n = 21 (I = 0, 1, 2, . . .) 

/ fv\ f A cosk n Y \Y\ <a 

and for odd n = 21 + 1 (I = 0, 1,2,.. .) 

, , f Bsmk n Y \Y\ < a , . 

^n{Y) = \ , x »-i . A-6) 

^ nV ; \ B(-l)— \Y\>a V ' 
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